Transconductance of a double quantum dot system in the Kondo regime 
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We consider a lateral double-dot system in the Coulomb blockade regime with a single spin- 1/2 on 
each dot, mutually coupled by an anti-ferromagnetic exchange interaction. Each of the two dots is 
contacted by two leads. We demonstrate that the voltage across one of the dots will have a profound 
influence on the current passing through the other dot. Using Poor Man's scaling, we find that the 
Kondo-effect can lead to a strong enhancement of this transconductance. 
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Nano-structured electron systems offer the unique pos- 
sibility to prepare and probe highly nonequilibrium sta- 
tionary states. Energy dissipation and therefore heat- 
ing, characteristic of any nonequilibrium state, takes 
place in the bulk reservoirs connected to the nano- 
structure, where it is absorbed causing a negligible 
change of state. The nature of such a nonequilibrium 
state may be conveniently probed by transport mea- 
surements, rather than thermodynamic or optical exper- 
iments (cf. e.g. Refs.[l|, Q)- In an ordinary quantum dot 
set-up the probing field, e.g. the applied bias-voltage, is 
also responsible for driving the system out of equilibrium. 
However, by separating the driving field and the probing 
field, many more possibilities for both studying nonequi- 
librium states and exploiting their properties arise. In the 
simplest realization of such a system one must have two 
parts of the system strongly coupled by some interaction 
and each contacted by two leads. Applying a sufficiently 
high bias across one part will then influence the current 
through the other part. This coupling is mediated by 
the change in nonequilibrium occupation numbers on the 
nano-structure caused by the driving field, which in turn 
governs the current-response to the probing field. 

In this letter we investigate the properties of a system 
of this type (cf. Fig. [T]), that is two quantum dots in the 
Coulomb blockade regime with odd electron occupation, 
giving rise to spins Sl,b. = 1/2, coupled by an exchange 
interaction K. Each dot is contacted by two leads, mod- 
eled by exchange tunneling amplitudes Jl,r ~ /Ec, 
in terms of lead-dot tunneling amplitudes t and charging 
energy Eq- This setup is very similar to the double-dot 
devices studied for example in Refs. where the 

central region was either a quantum-dot, a quantum-wire 
or a tunneling barrier, which could then support an ex- 
change interaction between the two spins. The following 
discussion docs not rely on the physical nature of the ex- 
change coupling. In practice, it can be due to a simple 
superexchange mechanism '6'| or an RKKY-intcraction, 



FIG. 1: (color online) Double Quantum Dot System: Two 
Kondo impurities denoted L, R each connected to two leads 
1, 2 and 3, 4, respectively, and coupled to each other by a spin 
exchange interaction K. 



as is most likely the case in Refs. [s, 4*], or even due 
to the so-called Kondo-doublct interaction suggested re- 
cently 01 . 

The theory of the two-inipurity Kondo model (2IKM) 
per se is well established |8|, but with the experimen- 
tal progress in quantum-dot systems, issues like the 
quantum critical properties remain an active field of 
research [l^, 11 1. The influence of the exchange- 
interaction K and an applied magnetic field on the non- 
linear conductance diL/dVL was studied in Refs. 12. l^. 
and it was demonstrated how the nature of the ground- 
state {Kondo-screened or local singlet) can be extracted 
from the temperature dependence of the linear conduc- 
tance. 

In this present setup with two sets of biased leads, we 
are facing an entirely new problem in which correlations 
and nonequilibrium effects must be treated side by side. 
We calculate the transconductance diL / dVu and demon- 
strate a transistor-effect in which a current through the 
left dot, II, as a response to a bias on the left dot, Vl, 
is switched on only when the bias over the right dot, Vji, 
exceeds a threshold given hy K > 0. This threshold ef- 
fect rests on the fact that the exchange-tunneling current 
from lead 1 to lead 2, 7^, is zero as long as the spins are 
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locked in a singlet. A sufficiently large bias across the 
right dot will, however, lead to a partial occupation of 
the triplet state, which does allow for a finite current 1^. 
As we shall demonstrate, the nonequilibrium polariza- 
tion due to Vr, and thereby the transconductance peak 
at Vfl ~ K, will be enhanced by the Kondo-effect. In the 
case of ferromagnetic coupling, the triplet ground state 
gives rise to a spin-1 Kondo-effect 12, ll3[. In this case, 
a bias-voltage over the right quantum dot will influence 
the corresponding zero-bias conductance peak and sin- 
glet sateUite peaks in diL/dVL, partly by decoherence 
and partly by a redistribution of the singlet/triplet oc- 
cupations. However, these effects are less dramatic and 
we therefore restrict our attention to the antiferromag- 
netic case, K > 0, which exhibits the afore mentioned 
threshold effect. 

The model. We model the double dot by two spin i 
operators Sl, Sij mutually coupled by exchange interac- 
tion K and each coupled to a pair of leads (couplings 
Jl,r > 0) as sketched in Fig. [T] The Hamiltonian reads 



(1) 



n,m— 1,2 n,m— 3,4 

where n, m = 1, 2, 3, 4 labels the leads, assumed to be in 
equilibrium at chemical potentials cl^^^^ creates elec- 
trons in lead n of momentum k and spin cr and the spin- 
operator s„m = J2k,k',cr,a' cik<T(^'J,'^'/2)cmk'<T' dcuotes re- 
spectively the conduction electron spin (n = m) and the 
exchange tunneling operator (n 7^ m). We do not allow 
for charge transfer between the two dots. Experimen- 
tally the charge transfer may be suppressed by a suitable 
arrangement of tunneling barriers (cf . Q ) . 

For K = 0, the two spins will undergo separate screen- 
ing by the conduction electron spins in each their set 
of leads as the temperature is lowered beneath the two, 
generally different, Kondo-temperatures Tk.l and Tk,r- 
However, at sufficiently large bias voltages cVl ^ Tk,l 
and/or cVr ^ Tk,r where eVL,eVff are the chemical 
potential differences of leads n = 1,2 or leads n = 3,4, 
resp ectively, these Kondo effects are partially suppressed 
[3[l^. A finite anti- ferromagnetic exchange coupling, 
K, on the other hand, will correlate the spins 8^,8^, 
and depending on the relative magnitude of K and the 
Kondo-temperatures, the two spins will either lock into 
a singlet or be screened by their respective lead elec- 
trons 0. We will be interested in the regime K ^ 
max{TR-,L, Tk,r}, when the Kondo effect cannot develop 
fully. In this case it is necessary to treat the exchange 
coupling K exactly, whereas it is possible to treat the 
Kondo effect perturbatively. 

The spin-states of the isolated double-dot are the sin- 
glet {S — 0) and triplet [S — 1) states of the two cou- 
pled spin 1/2. They form the appropriate basis states 
in this regime even when the coupling to the leads is 



turned on. For convenience, we employ a pseudo-boson 
(bond- operator [l6| ) representation of the four impurity- 

7 G {s,+,0,-} 
states are lUs — 
therefore K Sl 
constraint Q ; 



states, with creation operators {bl^} — {s\tl^,tl,t^_), 

The energy eigenvalues of the four 
-fx, Wo = L0+,- = = LOt, and 



u}^b\bj, 



4-- - wt, 
subjected to the 



+ V_t^ = 1. The con- 
straint may be exactly imposed by adding a term XQ 
to the Hamiltonian and taking the limit A ^ 00 17 1. 



In terms of the pseudo-bosons, the spin-1/2 operators 
are given by 5'£ ^ = ^{±sHo ± tjs + — t'Lt-) and 

sIr = ^,{sIr + ^si^) = {SZ,rV = U±^^t_T 

t^s + t^to + tji-) ^16], or in compact notation Sl.r = 
h '^■y 1' b\T^L,R;-if'bj' , whcrc Tl^h-j^i is a vector of three 
4x4 matrices, T^,Ty,T^. Notice that this exchange- 
tunneling model contains no operators of the form s^s, 
which immediately implies that the current 1^ is zero un- 
less transitions to the triplet-states become energetically 
allowed. 

It is known that away from the particle-hole symmetric 
point a finite potential scattering term will give rise to 
a current below the threshold. Nevertheless, we shall 
assume that a set of individual gates on each dot can 
be adjusted so as to make the system (nearly) particle- 
hole symmetric and hence neglect potential scattering 
altogether. In the setup by Sasaki et al. such tuning 
was made possible by two individual plunger-gates. 

Nonequilibrium perturbation theory. In equilibrium, 
i.e. for vanishing bias voltages, the states I7) of the dou- 
ble dot are thermally occupied, with occupation num- 
bers = {b\b^) = exp(-tj-y/T)/X]^' exp(-w^//r). In 
a stationary nonequilibrium situation with sufficiently 
large currents through the nano-structure, the occupa- 
tion numbers are determined by the currents, i.e. the 
couplings and the voltages, rather than temperature [l§ |. 
A further complication is that nondiagonal parts of the 
density matrix n^yi = {bl^bj') may be nonzero due to the 
mixing of states by the coupling to the leads. Confining 
our considerations to zero magnetic field this will not be 
the case and furthermore the triplet-states will be equally 
occupied, i.e. = nt, for 7 = -|-,0, — . It is useful to 
define the polarization p — ~ nt, and we note that 
(Sl ■ Sr) = -fp. 

The nonequilibrium occupation-numbers are deter- 
mined from the steady-state quantum Boltzmann equa- 
tion, '^^1^^ [Tj^inY — Ty-yn^] — 0, with transition-rates 
given by the golden rule expression 

TT °° 

^77' ^J2^ T^'7'T77' J (^)37'7 ('^-^7'7)^77' (^)' 

'^n,m —00 

(2) 

where {n,m) takes the values (1,2) or (3,4) and T"™ 
Ti or Th, uj^y - LO^-LJy and Fl^^pito) = /(w-a*„)[1- 
f{uj — fim — '^77')], fi'^) being the Fermi-function. The 
dimensionless exchange-couplings, g"™ = vfJ""^, (as- 
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suming same density of states, vp, in all leads) will only 
acquire frequency-dependence under renormalization as 
we shall discuss later. The Boltzmann equation is read- 
ily solved together with the constraint '^^n^ — 1, and 
using the three- fold degeneracy of the triplet, one finds 

thatp = (r,t-rt,)/(r,t-f3rt,). 

The current through the left dot is obtained as 



2 

iLiV) =Y JE"7TL;77' ' Tl;^^ (3) 
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and a corresponding expression holds for Iii{V), with 
leads 1, 2 replaced by 3, 4. 

From Eq. ([3]) it is clear that Vr only enters via the 
occupation numbers, and using the unrenormalized cou- 
plings the transconductance is given by 



diL 
dVp 



[F_(Vi) - cVl] ^F+{Vn), (4) 



where F±(F) = {K-eV)nB{K-eV)±{K + eV)nBiK + 
eV), in terms of the Bose-function ns(e) — l/(e'^/-^ — 1), 
and 3? = i9''9'^?/Emj9""'?- dlL/dVn is an odd 
function of Vr which diverges like T / {K—eVR) for cVr 
K. As we shall see below, this divergence is removed 
once the broadening of the excited triplet-states is taken 
into account, and one is left with a pronounced peak in 
the transconductance. Assuming that g™" — gL(R) for 
m,n & {1,2}, {3, 4}, one has g"^ = -I- gR^)~^ which 
shows that the transconductance peak is expected to be 
largest in a symmetric setup where g^ = gR- 

Renormalized perturbation theory. Since the higher or- 
der terms in the perturbation series are plagued by log- 
arithmic singularities characteristic of the Kondo effect, 
we apply Poor man's scaling to sum up the leading log 
arithms in this series. 



As demonstrated earlier 
this is conveniently done by keeping track of the gener- 
ated frequency dependence of the couplings. Within this 
approach the renormalized coupling functions are cal- 
culated from the following set of renormalization group 
(RG) equations describing the flow of couplings under a 
reduction of the bandwidth D: 



mg Ic 

Till 



dluD 



^ I 

+g':l(^ll+K)g\^{^ll)Q^^^,^K 
+ 2g\T{^ll)9u{m)Q^-^.^ 



(5) 



- -T.{9'^{M{pi)Q^^,^ 

+ g\i{pi~K)g[T{pii)Qu-^.,+K}, (6) 

with g'^i^iuj) = g(™"(cj-if). Here 6^ = Q{D^ -uj^ -T^) 
is the step function and F is a Korringa-like decoherence- 
rate due to particle-hole excitations in the leads. In gen- 
eral, F involves both self-energy, and vertex corrections 



FIG. 2: (color online) Polarization p = Us — nt versus voltage 
eVn/K for different values of K/Tk at T /Tk =0.1 and volt- 
age Vl = Inset shows the derivative of the curves in the 
main panel. 



and may depend on the intermediate state for which it 
describes the broadening (cf. "^j. Since, however, F 
appears only under the logarithm we will neglect these 
details here and simply use the maximum F = max{F^-y'} 
with F^^/ defined in Eq. ^ . The renormalized couplings 
and F are calculated self-consistently, then we solve for 
the occupation numbers and calculate the current us- 
ing the renormalized coupling- functions in Eq. ([3]). The 
perturbative RG approach is valid for K ^ T^, where 
Tx = ZJoe^^/^^^^"-* in terms of bare coupling g{DQ) and 
the half-bandwidth Dq. 

Renormalized conductances. In Fig. [2] we plot the po- 
larization p as a function of cVr/K at Vl = 0, for several 
values of K/Tk and at temperature T/Tk = 0.1. As cVr 
approaches K the triplet states become partially popu- 
lated and eventually p tends to zero as 1 /Vr. Since we use 
T <gC K, temperature is practically zero and the broad- 
ening is governed by K and Vr. In the inset of Fig. [2] 
we plot the derivative —Kdp/dVR vs. cVr/K, showing 
a pronounced peak which is increasingly smeared as the 
ratio K/Tk is made smaller and the Kondo-correlations 
become more pronounced. 

From the polarization p and the renormalized cou- 
plings, we can now determine the current from Eq. p|). 
As mentioned above, electron transport via exchange- 
tunneling always involves excitation out of the singlet 
ground state. This requires an energy of the order 
of K and therefore leads to a threshold behavior near 
cVl ^ K, as sketched also in Ref. [l2[. Alternatively, 
the singlet-triplet gap can be overcome by a finite Vr, 
which will lead to a finite dl^/dVL even for cVl < K. 
This is clearly seen in Fig. [3] When becomes larger 
than K, a zero-bias peak is seen to develop, the magni- 
tude of which is determined by the ratio TK/T^. In the 
limit of zero probing-voltage, we find that 



diL 

li™ ITT 

Vl^o dVL 



h [ilniVTK /TkW 



(7) 
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FIG. 3: (color online) Diflterential conductance cHl/cIVl in 
units of e'^/h for fixed T/Tk = 0.1 and l/ln{K/TK) = 0.2 
versus bVl/K at different values of cVr/K. Inset: dlh/dVL 
at 1/l — > versus cVr/K, full circles: perturbative RG, empty 




FIG. 4: (color online) Transconductance cHl/iIVr in units of 
e'^/h multiplied by (isT/eVi) versus eVn/K for fixed T/Tk = 
0.1, probing voltage eVl/K — 0.05 and different values of 
K/Tk- Inset: Transconductance for 1/ \n{K/TK) = 0.2 and 
different cVl/K. 



where both nt and F depend on the driving- vohage Vr. 
The inset of Fig. [3] shows the increase of the linear con- 
ductance with Vji, as calculated respectively by pertur- 
bative RG and 2"'^ order perturbation theory (PT). The 
initial rise with Vr is set by the triplet-occupation num- 
ber nt — {I — p)/4. The amplitude is logarithmically en- 
hanced due to the Kondo effect and falls off very slowly 
(F oc Vr) for Vr > K. 

Fig. 0] shows the transconductance versus cVr/K for 
Vl 0. Enhancing the Kondo-correlations by lowering 
K /Tk is seen to amplify the transconductance but also to 
smear the threshold. By scaling with K/ cVl we eliminate 
the 1/K dependence seen in Eq. ^ as well as most of the 
dependence on the probing field for cVl < K (cf. inset 
of Fig.SD. 

Conclusion. We studied the mutual influence of spa- 
tially separated currents flowing through two quantum- 
dots coupled only by spin-exchange. The strong interrela- 



tion between the currents and nonequilibrium occupation 
numbers for the singlet and triplet spin-states correlates 
the two currents, and we showed that this can be ob- 
served as a marked peak in the transconductance. Since 
the spin exchange interaction requires well-defined spins 
on the quantum dots the coupling to the leads should 
be relatively small, but opening up towards the Kondo- 
regime our perturbative RG calculations demonstrate 
that weak Kondo-correlations {Tk <C K) will in fact en- 
hance the transconductance peak. This pronounced sig- 
nal in the transconductance provides a new experimen- 
tal means of probing the exchange-coupling between two 
quantum-dot spins. Moreover, as evident from Eq. ([7]), 
this mode of operation provides direct experimental ev- 
idence for the voltage dependence of the nonequilibrium 
distribution fimction nt- 
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